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We derive hydrodynamic-like equations that are applicable to short-time scale color phenomena 
in the quark-gluon plasma. The equations are solved in the linear response approximation, and 
the gluon polarization tensor is derived. As an application, we study the collective modes in a 
two-stream system and find plasma instabilities when the fluid velocity is larger than the speed 
of sound in the plasma. The chromo-hydrodynamic approach, discussed here in detail, should be 
considered as simpler over other approaches and well-designed for numerical studies of the dynamics 
of an unstable quark-gluon plasma. 



I. INTRODUCTION 

Bulk features of electromagnetic plasmas are usually studied by means of fluid equations [1]. To get more detailed 
information, one refers to kinetic theory [1]. Since the fluid equations are noticeably simpler than the kinetic ones, 
the hydrodynamic approach is also frequently used in numerical simulations of the plasma evolution, studies of the 
nonlinear dynamics, etc. [1]. The situation in the field of quark-gluon plasma studies is rather different. Although 
chromo-hydrodynamics equations were discussed by several authors over a long period of time [2-12], they were not 
carefully studied. The field of their applicability was not established and very few important results were obtained 
by means of them. Consequently, the chromo-hydrodynamics has not attracted much attention. Instead, field theory 
diagrammatic methods have been successfully applied to reveal equilibrium properties of the quark-gluon plasma [13], 
while transport theory [14-16] proved efficient in the non-equilibrium studies. In particular, the important role of 
instabilities in the quark-gluon plasma evolution was clarified within the kinetic theory approach, see the review [17]. 
The kinetic equations were also a basis of extensive numerical simulations of the unstable QCD plasma [18-24]. The 
early stage of relativistic heavy-ion collisions, when the quark-gluon system is produced, was effectively studied using 
methods of classical field theory, see e.g. the review [25] and very recent publications [26,27]. 

Inspired by the success of hydrodynamic methods in the electromagnetic plasma, we discuss the approach to be 
applied to the quark-gluon plasma. Before going to the main subject of our study, however, a very important point has 
to be clarified. Real hydrodynamics deals with systems which are in local equilibrium, and thus it is only applicable 
at sufficiently long time scales. The continuity and the Euler or Navier-Stokes equations are supplemented by the 
equation of state to form a complete set of equations. The equations can be derived from kinetic theory, using the 
distribution function of local equilibrium, which by definition maximizes the entropy density, and thus the function 
cancels the collision terms of the transport equations. Such a real chromo-hydrodynamics was derived in [10] where the 
state of local equilibrium was found, using the collision term of the Waldman- Snider form. The chromo-hydrodynamics 
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has occurred trivial in the sense that ahhough the local equilibrium can be colorful, all color components of the plasma 
move with the same hydrodynamic velocity. Therefore, chromodynamic effects disappear entirely once the system is 
neutralized. It actually happens even before the local equilibrium is achieved [28]. Thus, there is no QCD analog of 
the magneto-hydrodynamics which is well known in the electromagnetic plasma. The magneto-hydrodynamic regime 
appears due to a large difference of electron and ion masses which effectively slows down a mutual equilibration 
of electrons and ions. Therefore, at a relatively long time scale one deals with the charged electron fluid in local 
equilibrium which moves in a passive background of positive ions. 

Since the hydrodynamic equations express the macroscopic conservation laws, the equations hold not only for 
systems in local equilibrium but for systems out of equilibrium as well. In particular, the equations can be applied at 
time scales significantly shorter than that of local equilibration. At such a short time scale the collision terms of the 
transport equations can be neglected. However, extra assumptions are then needed to close the set of equations, as 
the (equilibrium) equation of state cannot be used. In the electromagnetic plasma physics, several methods to close 
the system of equations were worked out. 

In this paper we discuss not a real hydrodynamics describing the quark-gluon plasma in local equilibrium, but the 
fluid equations which are valid at much shorter time scale. The approach is designed to study temporal evolution 
of the unstable QCD plasma. In Sec. Ill the fluid equations are derived from the kinetic theory which is presented 
in Sec. II. The equations are solved in the linear response approximation in Sec. IV and the polarization tensor is 
derived. As an example, the collective plasma modes in the two-stream system are analyzed in Sec. V. The paper is 
closed with a few remarks on the hydrodynamic approach. 

Throughout the paper we use the natural units with c = h = ks = ^ and the metric (1,-1,— 1,-1). 



II. KINETIC THEORY 



In this section we briefly present the transport theory of quarks and gluons [14,15] which is used to derive the fluid 

equations. 

The distribution function of quarks Q{p, x), where p is the four-momentum and x is the four-position, is a hermitian 
Nc X Nc matrix in color space (for a SU(Afc) color group). The distribution function is gauge dependent and it 
transforms under a local gauge transformation U as 

Q{p,x)^U{x)Q{p,x)UHx) . (1) 

The distribution function of antiquarks, which we denote by Q{p, x), is also a hermitian Nc x Nc matrix in color space 
it also transforms according to Eq. (1). The distribution function of gluons is a hermitian {N^ — 1) x (-^c ~ l) matrix, 
which transforms as 

G{p, x) U{x) G{p, x) U\x) , (2) 

where 

Uab{x) = 2Ty[t''U{x)t^U\x)] , (3) 

with r", = 1,..., Nc — 1 being the SU(7Vc) group generators in the fundamental representation with Tr (r^r'') = \5°'^. 

The four-momentum p, which is the argument of the distribution functions, is not constrained by the mass-shell 
condition but it is taken into account in the momentum integration measure. We use the notation 

/■-/l^^ewfe,^), (4) 

and the color current is expressed in the fundamental representation as 

fix) =-\j^P^ [Q^^^) - - ^Tr[Q(p,ar) - Q(p,ar)] + 2r«TV[T«G(p,a;)]] , (5) 

where g is the QCD coupling constant. A sum over helicities, two per particle, and over quark flavors Nj is understood 
in Eq. (5), even though it is not explicitly written down. The SU(A^c) generators in the adjoint representation are 
expressed through the structure constants T^" = —ifabc, and are normalized as Tr[T''T^] = Nc6'^^. The current can 
be decomposed as j^{x) = jli{x)T"- with j^{x) = 2Tr{Taj^{x)). 

The distribution functions of quarks, antiquarks and gluons satisfy the transport equations: 
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P^D^Q{p,x) + ^-p^ {F^.{x),d;Q{p,x)} = C[Q,Q,Gf\ , (6) 
p'^D^Q{p,x) - ^p^^ {F^,{x),d;Q{p,x)} =C[Q,Q,G], (7) 

P^V^Gip, x) + ^p^ {J'^Ax), d;G{p, x)} = Cg[Q, Q, G] , (8) 

where {...,...} denotes the anticommutator and dp the four- momentum derivative; the covariant derivatives and 
■D^ act as 

Dij,=d^- ig[Aij,{x), ... ] , = 9^ - ig[Aij.{x), ... ] , 

with and being four-potentials in the fundamental and adjoint representations, respectively: 

A^ix) = A^ix)T\ A'^ix)=T'^A'iix) . 

The strength tensor in the fundamental representation is F^,^ = d^A^ — dj^A^ — ig[A^, A^,], while J^^i, denotes the 
field strength tensor in the adjoint representation. C, C and Cg represent the collision terms which are neglected in 

our further considerations, as we are interested in short-timc-scale phenomena which are 'faster' than the collisions. 
The transport equations are supplemented by the Yang-Mills equation describing generation of the gauge field 

D^F^^^{x)=f{x), (9) 

where the color current is given by Eq. (5). 



III. FLUID EQUATIONS 

We assume here that there arc several streams in the plasma system, and then the distribution function of, say, 
quarks Q{p,x) can be decomposed as Q{p,x) = J2aQo:{p,x) where the index a labels the streams. Since we are 
interested in the short-time scale phenomena - 'faster' than collisions - the complete quark distribution function 
Q{p,x) as well as the distribution fimction of every stream Qa{p,x) satisfy the collisionless transport equation. The 
streams interact with each other only via the mean-filed which is generated in a self-consistent way by the current 
where every stream contributes. We note that the decomposition of the distribution function Q{p,x) into streams 
is, in principle, not unique. However, if the stream distribution functions Qa{p,x) solve the transport equations 
combined with the field generation equation, the complete distribution function solves the transport equation and the 
field generation equation as well. 

Further analysis is limited to quarks but inclusion of anti-quarks and gluons is straightforward. Integrating the 
collisionless transport equation (6) satisfied by Qa over the four- momentum with the integration measure (4), one 
finds the covariant continuity equation 

= , (10) 

where n'^ is A^c x -^c matrix defined as 

K{x)= ( p>'Qaip,x). (11) 
Jp 

The four-flow transforms under gauge transformations as the quark distribution function, i.e. according to Eq. (1). 

Multiplying the transport equation (6) by the four- momentum and integrating the product with the measure (4), 
we get 

I)^Tr-f{F/,<} = 0, (12) 

where the energy-momentum tensor is 

Tii^{x)^ I p^^p^Q^{p,x) . (13) 
Jp 

We further postulate that and Tj^" have the form of an ideal fluid i.e. 
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TH^ix) = l{e„{x) +p4x)){K{x),K{x)} -Pa{x)g' 



(14) 



(15) 



where the hydrodynamic velocity is, as Ua, and pa, sl Nc x Nc matrix. The anticommutator of ttJJ and is 
present in Eq. (15) to guarantee the symmetry of T^" with respect to /i ^ z/ which is evident in the definition (13). 

Since Eqs. (14, 15) are of crucial importance for our further considerations, let us discuss them in more detail. 
There are two aspects of the ideal-fluid form of n(;J and T^"^ to be justified: the Lorentz structure and the matrix 
structure. We start the discussion with the Lorentz structure. Eqs. (14, 15) assume that and T^'^ can be expressed 
as products of the Lorentz scalars tIq,, Cq, and Pa and the four- vector m^. Such a Lorentz structure appears when 
the distribution function Qa{p,x), which gives and T^"^ through Eqs. (11, 13), is isotropic in the local rest frame 
where = (1,0,0,0). One does not need to assume that the fimction is of local equilibrium form, the isotropy is 
a sufficient condition. And we note that the isotropy of the momentum distribution, required Eqs. (14, 15), is not a 
serious limitation because we actually consider not a single stream but a superposition of several streams. And then, 
one can easily model an arbitrary momentum distribution as a sum of several momentum distributions, everyone 
isotropic in its Lorentz rest frame. Thus, the Lorentz structure expressed by Eqs. (14, 15) seem to be well justified. 

The problem of the matrix structure of Eqs. (14, 15) is more complicated. We first note that and T^^ have to 
be matrix value functions as the distribution function Qa{p,x) which defines them through Eqs. (11, 13) is a matrix 
in a color space. The question is whether n„, £„, Pa and u'^ have to be all matrix quantities. When the system under 
consideration is colorless, these quantities can be treated as scalars in the color space. However, as explicitly shown 
in Sec. IV, even small perturbations of these quantities induced by a color field are of matrix value. Therefore, Uq,, 
Ca, Pa and have to be all treated as matrices. 

When and T^^ are expressed through n^, £„, Pa and u'^, one faces the problem of ordering the matrices. One 
can define in three ways: = no-u^, = and = {ua, 12 where {...,...} is the anticommutator. 

The three diS'erent definitions obviously lead to three different fluid equations coming from D^n^ = 0. We first note 
that the problem of ordering does not appear when small perturbations of colorless state are studied because the 
linearized equations are all the same, as follows from considerations presented in Sec. IV. The problem of ordering is 
also absent for systems, which are significantly colorful, if the matrices belong to the Cartan subalgebra of SU(7Vc), 
as then they commute with each other. Unfortunately, we have not found a general solution of the ordering problem. 
Numerical simulations will be very helpful to clarify whether the three formulations are equivalent to each other. 

In the case of an Abclian plasma, the relativistic version of the Euler equation is obtained from Eq. (12) by removing 
from it, following Landau and Lifshitz [29], the part which is parallel to u^. An analogous procedure is not possible for 
the non-Abelian plasma because the matrices n„, w^, and u"^^ in general, do not commute with each other. Thus, one 
has to work directly with Eqs. (10, 12) with and T^^ of the form (14, 15). The equations have to be supplemented 
by the Yang-Mills equation (9) with the color current of the form 



where only the quark contribution is taken into account. 

The fluid equations (10, 12) do not form a closed set of equations even when the chromodynamic field is treated 
as an external one. There are five matrix equations {v = 0, 1,2,3) and six unknown matrix functions: na{x), ea{x), 
Pa{x) and three components of the four- velocity u^{x). We note that the constraint u^(x)wq,^(x) = 1 is imposed. 

There arc several ways to close the system. In the case of real hydrodynamics, which describes the system in 
local thermodynamic equilibrium, one adds the equation of state. Although the system under consideration is not in 
equilibrium, we can still add a relation analogous to the equation of state. The point is that the energy-momentum 
tensor (13) of a weakly interacting gas of massless quarks is traceless {T^f^ix) = 0) because = 0. Then, Eq. (15) 
provides the desired relation 



due to the constraint u'^{x)ua nix) = 1. In fact, the complete energy-momentum tensor of any conformal theory 
is traceless but quantum effects, as the trace anomaly, can spoil conformal invariance. The phenomenon of mass 
generation in a medium modifies Eq. (17) but the effect is small in the pcrturbativc regime studied here. We also 
note that {x) and Pa {x) are not the energy density and pressure of the equilibrium system but the matrices in color 
space which in equilibrium become the energy density and pressure, respectively. 

Another method to close the system of equations is to neglect the gradients of pressure (pq,(.t)), assuming that the 
system's dynamics is dominated by a self-consistently generated chromodynamic field. In the following section, where 
the equations are solved in the linear response approximation, we use both methods to close the system. 




a 



(16) 



eaix) = 3pa{x) 



(17) 
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Our analysis is limited to quarks but, as already mentioned, inclusion of antiquarks and gluons is straightforward. 
Since the distribution functions of quarks, antiqiiarks and gluons of every stream are assiimcd to obey the coUisionless 
transport equation, we have a separated set of fluid equations for quarks, antiquarks and gluons of every stream. 
The equations are coupled only through the chromodynamic mean field. The quarks belong to the fundamental 
representation of the SU(iVc) group and thus, the hydrodynamic quantities na{x), €a{x), Pa{x) and u'^{x) are Nc x Nc 
matrices. Antiquarks can be treated as in kinetic theory that is as belonging to the fundamental representation even 
so they belong, strictly speaking, to the transposed fundamental representation. Since gluons belong to the adjoint 
representation, the hydrodynamic quantities, and consequently the fluid equations, are {N^ — 1) x (-^c ~ -'-) matrices. 

IV. LINEAR RESPONSE ANALYSIS 

In this section the hydrodynamic equations (10, 12) are linearized around an stationary, homogeneous and colorless 
state described by n, e, p and u^^. We mostly skip here the index a to simplify the notation. The index is restored in 
the very final formulas. 

Because every stream is assumed to be colorless, the matrices n, e, p and are proportional to the unit matrix in 
the color space. Thus, 

fiu''- -^Ti[nu^] = , (18) 

which means that the color four-current of every stream vanishes. The quantities of interest are decomposed as 

n{x) = fi + Sn{x) , e{x) = e + Se{x) , (19) 

p{x) =p + 6p{x) , ut^ix) =u^' + Suffix) . (20) 

Since the state described by n, e, p and is assumed to be stationary, homogeneous and colorless, we have 

D''n = 0, D''e = 0, L>^p = , = . (21) 

And because we consider only small deviations from the stationary, homogeneous and colorless state, the following 
conditions are obeyed 

n > (5n , e > (5e , p > (5p , > 6u^ . (22) 

Actually, 6n, 6e, 6p and 6u^ should be diagonalized to be comparable to the n, e, p and m^. 
Substituting the linearized and T^", which are 

n^" = nu^' + nSu^' + Snu'' , (23) 

T^"' = (e + p)u^' u" -p g^"' + {Se + 5p)u^' u'' + {e + p){u^' SvT + Ju" u") - Sp g^"" , (24) 
into Eqs. (10, 12), one finds 

n Du^Sui^ + {D^5n) = , (25) 

{Df,{6€ + Sp))u^' u" + {e+p){u^' {D^Su") + {D^Su^") u") - D''5p - gF>"'nUi, = . (26) 
Projecting Eq. (26) on W^, one finds 

u'^D^Se + (e + p)D^6u'' = . (27) 
To derive Eq. (27) one has to observe that u^Ufj, = 1 = u^u^, and consequently, u^Su^ = 0[{6u)^) and u^D^dun = 

Acting on Eq. (26) with the projection operator {g^^ — u^Uv), one gets the linearized relativistic Euler equation 

(e + p)u^Di'5u'' - {D" - u''u^D'')Sp - gnUf.F'"' = . (28) 

As already mentioned, Eqs. (25,27,28) do not form a closed set of equations. In the following two sections we solve 
the fluid equations (25,27,28), adopting two different methods to close the system. 
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A. Pressure gradients neglected 



The system is closed when the pressure gradients are neglected, which physically means that the system's dynamics 
is dominated by the mean field. When the term with the pressure gradient is dropped in Eq. (28), the equation 
(27) effectively decouples from the remaining equations, and one has to solve only two equations (25, 28) where 6e is 
absent. 

In principle, Eqs. (25, 28) with the pressure term neglected can be formally solved in a gauge covariant manner, using 
the inverse operator of the covariant derivative. However, we arc interested here only in computing the polarization 
tensor. Thus, we replace the covariant derivatives by the normal ones to fully linearize the equations, and the gauge 
independence of the result is checked a posteriori. After performing the Fourier transformation, which is defined as 

J (04e-*v(fc), (29) 

Eqs. (25, 28) get the form 

u'^kf.Sn + nk^Su^ = , (30) 

«(e + p)u^k''Su'' + gnu^F'"' = . (31) 
Eqs. (30, 31) are easily solved providing 

Sn = -ig^_p^F^-^ , (32) 
e + p (u • k)^ 

= ig _^ F''" . (33) 
e + p u-k 

Since Tr[F^^] = 0, the induced current equals 

a 

where the index labelling the streams present in the plasma system is restored. Because the linearized strength tensor 
equals F>"'{k) = -ik^" A" {k) + ik"" A^" [k) , one finally finds 

53^ {k) = -W{k)A,{k) , (35) 

with 

n (fc)--Yl.|-^ ^-^fcp • (36) 

The polarization tensor 11'*'' (fc) is proportional to the unit matrix in color space, it is symmetric and transverse 

(k^Ii^'^{k) = 0), and thus it is gauge independent. 

The polarization tensor obtained within the kinetic theory in analogous approximation [30] is 

W^{k) = --JjiP) (P • ^)(^'^^'^ + ^'^^^) - ''^P'^P'^ - ■ , (37) 

where f{p) is the distribution function of quarks in the colorless, stationary and homogeneous state. One observes 
that Eq. (37) transforms into Eq. (36) when 

m = J2n^ul (5(3) (p _ h±P^ . (38) 

The delta-like distribution function (38), and consequently the approximation, where the pressure gradients are 
neglected, hold for systems where the thermal momentum (pthermai) is much smaller than the collective momentum 
(PcoUcc) of the hydrodynamic flow. In the electromagnetic plasma of electrons and ions, such a situation occurs for a 
sufficiently low temperatures when the effects of pressure are indeed expected to be small. For the system of massless 
partons, the condition pthermai ^ PcoUec is achieved by requiring that PcoUec is large rather than pthermai is small. For 
massless particles in local equilibrium pthcrmai ^ T, where T is the local temperature, while the formula (38) clearly 
shows that PcoUec ~ TjaVa where Va and 7q. are the velocity and Lorentz factor of the collective flow. Therefore, the 
condition PcoUec > Pthermai requires 7a » 1- 
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B. Effect of pressure gradients included 



The equation (17) allows one to close the system of fluid equations, not neglecting the pressure gradients. Using 
the relation (17), one has to solve three equations (25,27,28). 

Performing the linearization analogous to that from the previous section and the Fourier transformation, the 
equations (25,27,28) get the form 



u^k^5n + nk^5u^ = 



(39) 



(40) 



Substituting obtained from Eq. (40) into the Euler equation (41), one gets 



3(m • k) 



^ {k'^k'' -vrk^'{u-k))\5u^ = i^g 



4^e(n-fc) 



Observing that 



[kakv — Uakviu ■ k)) [k^k^ - u''k^{u ■ k)) cx {k^k^ - Uak^{u ■ k)) , 
the operator in the left- hand- side of Eq. (42) can be inverted as 



1 



fc2-h2(w-fc)2 

and Eq. (42) is exactly solved by 



{k(jki/ Ufjki/iu ' fc)^ 



+ 



1 



3{u-ky 



-{k^k^" -u^k^iu-k)) 



5Ua 



.3 

^5 



1 



4^ e(u- fc) L^''^ fc2 + 2(u-fc)2 
Substituting 5n and Ju^ into the current (34), one finds the polarization tensor as 

1 



4ea (Wa • kf 



(w„ • fc)(fc''u^ + fc^w^) - (w„ • kfg^''' - k'^ui'^ul 

{ua ■ k)k^{ki^u'^ + k'^u^) - (m„ • kfk^'k'' - k'^u'^u'^ 
k"^ + 2(u„ • fc)2 



(41) 



(42) 



(43) 



(44) 



(45) 



where the stream index a has been restored. The first term in Eq. (45) corresponds to the polarization tensor (36) 
found when the pressure gradients are neglected while the second term gives the eff'ect of the pressure gradients. The 
first term is, as already mentioned, symmetric and transverse. The second term is symmetric and transverse as well. 
Thus, the whole polarization tensor (45) is symmetric and transverse. 

As explained below Eq. (38), the pressure gradients can be neglected and the distribution function can be approx- 
imated by the delta-like form (38) when 7^ ^ 1. Because the four- velocity is expressed as = (jaija^a), it is 
easy to check that the second term in Eq. (45) is small, when compared to the first one, for 7q ^ 1. However, one 
has to be careful here: when the wave vector k is exactly parallel to Vq and both are, say, in the z direction, then 
the two terms of the non-zero components of the polarization tensor (45) (H^", H''^, H^^) vanish as I/72. For the 
remaining non-zero components of (45) (11^^, 11^^), the second term is identically zero. We return to this point when 
the collective modes are discussed. 
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V. COLLECTIVE MODES IN THE TWO-STREAM SYSTEM 



As an application of the polarization tensors derived in the previous section, we discuss here collective plasma modes 

in the two-stream system. Within kinetic theory such a system was studied in [31 34]. The two-stream configuration 
has not much to do with the experimental situation in heavy-ion collisions where the momentum distribution of the 
produced partons is very different from the two-peak shape characteristic for the two-stream system. Thus, it should 
be treated only as a toy model. However, the model dynamics is rather nontrivial due to several unstable modes. 
While the two-stream configuration is rather of academic interest, one can model experimentally relevant situations 
with several streams. And in the linear regime, such a multi-stream configuration can be studied along the lines 
presented below. 

Since Fourier transformed chromodynamic field A^{k) satisfies the equation of motion as 



the general dispersion equation of the collective plasma modes is 







det 



k^g"" -k^k"" -W^ik) 







(46) 



(47) 



Due to the transversality of W^'^{k) not all components of 11'^'' (fc) are independent from each other, and consequently 
the dispersion equation (47) , which involves a determinant of a 4 x 4 matrix, can be simplified to the determinant of a 
3x3 matrix. For this purpose one usually introduces the color dielectric tensor e*-'(fc). where the indices i,j = 1,2,3 
label three-vector and tensor components, which is expressed through the polarization tensor as 

e'^{k)=S'^ + \w^k) , 



where k= (a;,k). Then, the dispersion equation gets the form 



det 



-k'k' -LJ^s'^k) 



0. 



(48) 



The relationship between Eq. (47) and Eq. (48) is most easily seen in the Coulomb gauge when = and k- A(fc) = 0. 
Then, E = iujA and Eq. (46) is immediately transformed into an equation of motion of E(fc) which further provides 
the dispersion equation (48). 

The dielectric tensor given by the polarization tensor (36), which neglects the effect of pressure, is 

,2 , 



ea+Po 



vl¥ + vik' (w^ - l<i^)vlvi 



w — k • 



(W - k • Va)^ 



(49) 



where Vq is the hydrodynamic velocity related to the hydrodynamic four- velocity u^; u)p is the plasma frequency 
given as 



2 _ 9 

U}„ = — ^ 

The dielectric tensor given by the polarization tensor (45), which includes the effect of pressure, is 



(w-k. 



(50) 



(51) 



(w - k • v„)(a;2 - k2)(5;jJ + 5^ F) - (cj - k • v„)2fc«F - (w^ _ X^fv^vi 



(a; - k • v„)2(a;2 - k2 + 27^(0; - k • v„)2) 



with the plasma frequency given as 



3<? 



(52) 



In the following subsections we consider the collective modes in the two-stream system with the wave vectors 
perpendicular and parallel to the hydrodynamic velocity. The velocities v^^ are chosen to be oriented along the axis 
z, i.e. Vet = (0, 0,z;a)- The index a, which labels the streams has two values, a = ±. For simplicity we also assume 
that 



V = V+ = —V- 



n = n+ = n_ 



e = e+ = e_ 



P = P+=P- 



(53) 



Then, the plasma frequency equals Up = g^r? j (e -|- p) or = 2>g^r? j (4e). 



8 



A. k_L V 



The wave vector is chosen to be parallel to the axis a;, k = (fc,0,0). Due the conditions (53), the off-diagonal 
elements of the matrix in Eq. (48) vanish. 



1. No pressure gradient effect 
With the dielectric tensor given by Eq. (49), the dispersion equation is 

(^2 _ ^2)(^2 _ ^2 _ ^2)^^2 _ ^2 _ ^2 _ ^2 ^l^) = Q , (54) 

where = WpU^. As solutions of the equation, one finds the stable longitudinal mode with = and the stable 
transverse mode with = uj^ + k"^. There are also transverse modes with 

^± = ^('^p - A' + fc' ± ^{^l - A2 + fc2)2 +4A2fc2) . (55) 

As seen, > but < 0. Thus, the mode W-i- is stable and there are two modes with pure imaginary frequency 
corresponding to w?. < 0. The first mode is overdamped while the second one is the well-known unstable Weibel 
mode leading to the filamcntation instability. 

As discussed below Eq. (38), the pressure gradients can be neglected when 7 1, and thus, the solutions (55) are 
physically relevant in this limit. Therefore, we write down the solutions (55) for 7 l^- 1 which are 

= i(A:2±.^A;4+4a;2fc2j . (56) 
As seen in Eq. (56), the modes are independent of 7, if it is sufficiently large. 



2. Effect of pressure gradients included 

As in the 'pressureless' case, there are two stable modes oS^ = Up and the stable transverse mode with u"^ = uip + k"^. 
The transverse modes corresponding to those given by Eq. (55) are obtained by solving the equation 

- .\k^ + .1) + X\.- - k-) - >^ \,^f~l'l-,^, = . (57) 

Defining the dimensionless variables a = oj/wp and b = k/cop, the above equation can then be rewritten as follows 

[a^ - F{b, v)a'^ - G{b, v)] = , (58) 



where 



Solutions of Eq. (58) are given by a = and 

= ^{F ± ^/F^TAG) (60) 
The solution with describes the unstable mode (a?_ < 0) provided G is positive. The condition G > gives 

v^>- and k^ < ujI . (61) 



As seen, the instability appears when the stream velocity is larger than the speed of sound in the ideal gas of massless 
partons and when the wavelength is sufficiently long. We also note that in the limit 7 » 1 the solutions (60) reproduce, 
as expected, those given by Eq. (56). 
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B. k II V 



The velocities and the wave vector are chosen to be oriented along the 2— axis i.e. v± = (0, 0, ±v) and k = (0, 0, k). 
Then, the matrix in Eq. (48) is diagonal. As in the case k _L v, wc discuss separately the collective modes given 
by the dielectric tensors (49) and (51). However, it is physically rather unjustified because, as explained at the very 
end of Sec. IV, the tensor (49) cannot be treated as an approximation of (51) even for j ^ 1. The effects of the 
mean-field and of the pressure gradients appear to be equally important. Therefore, only the results of subsection 
V B 2 are physically reliable. 



1. No pressure gradient effect 



With the dielectric tensor given by Eq. (49), the dispersion equation reads 

There arc two transverse stable modes with = + fc^. The longitudinal modes are solutions of the equation which 
can be rewritten as 

where coq = ujp/2j'^. With the dimensionless quantities x = lj/ujq, y = kv/uQ, Eq. (63) is 

{x^ - y^f - 2x^ -2y^ = 0, (64) 

and it is solved by 

x%=y^ + l± V4y2 + 1 . (65) 

As seen, a;^ is always positive and thus, it gives two real (stable) modes; is negative for < y < -\/2 and then, 
there are two pure imaginary modes. The unstable one corresponds to the two-stream electrostatic instability. 



2. Effect of pressure gradients included 

The transverse modes are the same as in the 'pressureless' case while the longitudinal modes are solutions of the 
equation 



1 1 



° V (w - kvf (w + kvf 



i^p ( {ujv - fc)2 [ujv + kf 



2 V(w-fct^)2(272(a;-A;t;)2+w2-A;2) {u + kv)'^{2^^{uj + kv)"^ + to"^ - k"^) 
which, in terms of the dimensionless quantities introduced in the previous section, is 



(66) 



1 ^fi + i^^-vM' \ + {v^ + y/v? \ (67) 

{x-yr\ ^ 2{x-y)^ + ^{x^-yyv^)J {x + y^y ^ 2{x + y)^ + j,{x^ - y^v^) J ^ ' 

The equation can be converted to 

{x^ - y^f [x^ - B{y, v)x'' - C{y, v)] = (68) 

with 
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which is solved by = and 

x4 = i(i?±VB2 + 4C) . (70) 

As for the case when the wave vector is perpendicular to the velocity, there are only instabilities when the stream 
velocity is bigger than the speed of sound in the ideal gas of massless partons, as C > is satisfied when 

^^'>^ and fc^olz^a;^ (71) 

It is also interesting to note that for the particular case u^ = l/3, C = 0, and the solutions of the dispersion relations 
get a very simple form. Namely, we have w?. = and 

0,^ = ^(0,^ + ^^). (72) 



VI. CONCLUDING REMARKS 

The fluid equations (10,12) with and T^" given by the formulas (14,15) form a closed set of equations when 
supplemented by the relation (17). If the chromodynamic field is generated self-consistently by the color fluid, the 
Yang-Mills equation (9) with the color current of the form (16) must be added. The derivation of these equations is 
the main result of our paper. Although some information is lost when going from kinetic theory to hydrodynamics, 
the dynamical content of the fluid equations is still very rich. The results of Sec. V, where the two-stream system is 
studied, demonstrate it convincingly. The system is shown to be unstable with respect to the magnetic and electric 
modes, if the hydrodynamic velocity is sufficiently large. 

As noted in the Introduction, the hydrodynamic approach is frequently used in studies of electromagnetic plasmas. 
Numerical simulations of the fluid equations are much simpler than those of kinetic theory because the distribution 
function depends on four-position and momentum while in the fluid approach one deals with a few functions of four- 
position only. For this reason we expect that the equations derived here will be useful in numerical studies of the 
unstable quark-gluon plasma. Hopefully, the late stage of instability development, when the dynamics is strongly 
nonlinear, can be attacked within the fluid approach. In the case of the two-stream system, our analytic results from 
Sec. V can be used to gauge numerical calculations. While the two-stream configuration is rather academic, one 
can model situations, which are relevant for quark-gluon plasma from relativistic heavy-ion collisions, with several 
streams. And then, our fluid equations can be applied. 

We also note that the hydrodynamic approach presented here can be extended to include the effect of collisions which 
can be modeled with transport coefficients (viscosities, conductivities) and additional terms in the fluid equations. 
Then, the approach would be applicable to longer time scales. 

Our derivation of the chromo-hydrodynamic equations is based on the parton kinetic theory valid in the regime 
of very weak couplings or, equivalently, of very high temperatures T » Tc, where is the temperature of the 
deconfinement phase transition. One may wonder what happens in the regime of not so weak couplings. Because the 
hydrodynamic equations are the expressions of the conservation laws, one can hope that the structure of the equations 
survives when the plasma is no longer weakly coupled. Lattice studies indicate that an equilibrium quark-gluon plasma 
at temperatures close to is very different from the ideal gas of massless quasiparticles [35]. The deviations can 
be, at least partially, accounted for in a quasiparticle model when the gluon dispersion relation is modified, see e.g. 
[36-38]. In such a case, our approach can be easily adapted for by changing the relation (17). We hope to report 
about those results in the near future. 



VII. ACKNOWLEDGMENTS 

We are grateful to Massimo Mannarelli and Mike Strickland for fruitful discussions. CM. was supported by MEG 
under grants FPA2004-00996 and AYA 2005-08013-C03-02, and by GVA under grant GV05/164. 



11 



[1] N.A. Krall and A.W. Trivolpiccc, Principles of Plasma Physics (Mc-Graw-Hill, New York, 1973). 

[2] K. Kajantic and C. Montonen, Phys. Scripta 22, 555 (1980). 

[3] D. D. Holm and B. A. Kupershmidt, Phys. Lett. A 105, 225 (1984). 

[4] D. D. Holm and B. A. Kupershmidt, Phys. Rev. D 30, 2557 (1984). 

[5] St. Mrowczynski, Phys. Lett. B 202, 568 (1988). 

[6] St. Mrowczynski, in Quark-Gluon Plasma edited by R. Hwa, Adv. Ser. Direct. High Energy Phys. 6, 185 (1990). 

[7] J. Bhatt, P. Kaw and J. Parikh, Phys. Rev. D 39, 646 (1989). 

[8] R. Jackiw, V. P. Nair and S. Y. Pi, Phys. Rev. D 62, 085018 (2000). 

[9] B. Bistrovic, R. Jackiw, H. Li, V. P. Nair and S. Y. Pi, Phys. Rev. D 67, 025013 (2003). 
[10] C. Manuel and St. Mrowczynski, Phys. Rev. D 68, 094010 (2003). 
[11] J. Dai and V. P. Nair, arXiv:hep-ph/0605090. 

[12] B. A. Bambah, S. M. Mahajan and C. Mukku, arXiv:hep-th/0605204. 

[13] J. P. Blaizot and E. lancu, Phys. Rept. 359, 355 (2002). 

[14] H. T. Elze and U. W. Heinz, Phys. Rept. 183, 81 (1989). 

[15] St. Mrowczynski, Phys. Rev. D 39, 1940 (1989). 

[16] D. F. Litim and C. Manuel, Phys. Rept. 364, 451 (2002). 

[17] St. Mrowczynski, Acta Phys. Polon. B 37, 427 (2006). 

[18] A. Rebhan, P. Romatschke and M. Strickland, Phys. Rev. Lett. 94, 102303 (2005). 

[19] A. Dumitru and Y. Nara, Phys. Lett. B 621, 89 (2005). 

[20] P. Arnold, G. D. Moore and L. G. Yaffe, Phys. Rev. D 72, 054003 (2005). 

[21] A. Rebhan, P. Romatschke and M. Strickland, JHEP 0509, 041 (2005). 

[22] P. Arnold and G. D. Moore, Phys. Rev. D 73, 025006 (2006). 

[23] A. Dumitru, Y. Nara and M. Strickland, arXiv:hep-ph/0604149. 

[24] P. Romatschke and A. Rebhan, arXiv:hep-ph/0605064. 

[25] E. lancu and R. Venugopalan, in Quark-Gluon Plasma 3, edited by R.C. Hwa and X.N. Wang (World Scientific, Singapore, 
2004). 

[26] P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96, 062302 (2006). 

[27] T. Lappi and L. McLerran, Nucl. Phys. A 772, 200 (2006). 

[28] C. Manuel and St. Mrowczynski, Phys. Rev. D 70, 094019 (2004). 

[29] L. D. Landau and E.M. Lifshitz, Fluid Mechanics (Porgamon, Oxford, 1963). 

[30] St. Mrowczyriski and M. H. Thoma, Phys. Rev. D 62, 036011 (2000). 

[31] Y. E. Pokrovsky and A. V. Selikhov, JETP Lett. 47 (1988) 12 [Pisma Zh. Eksp. Teor. Fiz. 47 (1988) 11]. 
[32] Y. E. Pokrovsky and A. V. Selikhov, Sov. J. Nucl. Phys. 52, 146 (1990) [Yad. Fiz. 52, 229 (1990)]. 
[33] Y. E. Pokrovsky and A. V. Selikhov, Sov. J. Nucl. Phys. 52 (1990) 385 [Yad. Fiz. 52 (1990) 605]. 
[34] St. Mrowczynski, Phys. Lett. B 214, 587 (1988). 
[35] F. Karsch, Nucl. Phys. A 698, 199 (2002). 

[36] M. I. Gorenstein and S. N. Yang, Phys. Rev. D 52, 5206 (1995). 

[37] D. Zwanziger, Phys. Rev. Lett. 94, 182301 (2005). 

[38] M. Bluhm, B. Kampfer and G. SofF, Phys. Lett. B 620, 131 (2005). 



12 



